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Abstract 

We show that recent apphcations of the homotopy perturbation method the Ado- 
mian decomposition method and the variational iteration method are completely 
useless for the treatment of nonlinear problems. 



1 Introduction 

In a series of papers we have shown that some particular applications of some 
variational and perturbation approaches (VAPA) like, for example, the homo- 
topy perturbation method (HPM), homotopy analysis method (HAM), Ado- 
mian decomposition method (ADM) and variational iteration method (VIM) 
are utterly useless for the study of nonlinear systems, even for the simplest 
models [1-7]. Unfortunately some journals do not accept criticisms of the pa- 
pers they publish. 

The purpose of this paper is the analysis of some recent applications of VAPA 
to simple models of nonlinear phenomena published in this journal. In Sec. [2] 
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we discuss the application of HPM to an exactly solvable Riccati equation. In 
Sec. [3] we analyse the application of the same approach to an exactly-solvable 
model for the interaction between a prey and a predator. In Sec. IHwe discuss 
the application of other VAPA to the simplest epidemic model. Finally, in 
Sec. |5] we draw some conclusions on the achievements and future of those 
VAPA. 

2 Homotopy perturbation method for a quadratic differential equa- 
tion 

Lately, there has been great interest in the application of analytical approxi- 
mate methods to the solution of several problems of nonlinear dynamics [8-13], 
some of which lead to time-power series [8-11]. Some of the models consid- 
ered by those authors are suitable for the study of prey-predator interac- 
tions [8,9,11-13]. 

It is curious that most of those authors resort to a power series description of 
the dynamics of the system because it is known that such approach is unable 
to provide a reasonable overall picture of the evolution which is what really 
matters in the case of, for example, prey-predator models. A time-power 
series is limited to a neighbourhood of the origin (for example, initial species 
population) by its convergence radius that is determined by the singularity 
closest to that point. Besides, most equations for nonlinear dynamics generate 
singularities spontaneously that move when the initial conditions change. We 
have already pointed out the limitations of several of those approaches in 
earlier communications [1-7]. 

The purpose of this section is to discuss some of those features of the nonlinear 
dynamics by means of the simple Riccati equation solved by Abbasbandy [10] 
by means of the HPM. 
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The model chosen by Abbasbandy [10] 
dY{t) 



dt 



2Y{t) - Y{ty + 1, r(o) = 



for the apphcation of the HPM is useful for present discussion because the 
Riccati equation can be solved exactly: 



One clearly appreciates that there is a pole in the complex t-plane that limits 
the convergence of the time-power series to t < \tc\, where 

^ % (v^+ l)l' + 7r2 ^ 1.274 (3) 



In other words, the HPM proposed by Abbasbandy [10] will be useless for 
t > \tc\ disregarding the order of the perturbation approach. In particular, the 
resulting time-power series will not reveal the stationary point Yg = \/2+l [14] 
to which the solution approaches a.s t oo. Such stationary points are most 
important for any physical description of a dynamical system [15,16]. 

The solution to the Riccati equation ([T]) with an arbitrary initial condition 
Y{0) = Yo is 

Y(t) = ^ ^-^ ^ ^ (4) 

(Fo + V2 - l) e2v^* - + 72 + 1 



and the radius of convergence of the time-series is 

^In (^ ^0- v^- 1 
4 ''[yo + V2-1 



^lnr°-f-^l (5) 



Notice that the singularity closest to the origin of the complex t-plane moves 
as the initial condition changes as mentioned above. The HPM performs even 
more poorly if Yq > for example, if Iq = 5 then |tc| ~ 0.261. 



3 



Prom the discussion above we clearly appreciate that the HPM proposed by 
Abbasbandy [10] (and also the time-power series produced by any other ap- 
proach [8,9, 11]) is unable to reveal the main features of nonlinear dynamical 
models. 

Straightforward application of the ADM [8, 14] also produces a time-power 
series, however, El-Tawil et al [14] proposed a multistage ADM that basi- 
cally leads to the expansion oiY{t + At) in Ai-power series. The radius of 
convergence of this series is given by 



|Aic|(i) = ^y41n(V2 + l) -8V21n(V2 + l)i + 8i2 + 7r2 (6) 

> — ~ 1-11 
4 

when Yq = 0. El-Tawil et al [14] chose AT ^ 1.11 and their multistage 
ADM yielded accurate results of Y{t) even for values of t sufficiently large 
that Y{t) ~ Ys. However, this approach is more closely related to numerical 
integration algorithms such as Runge-Kutta than to the analytical methods 
discussed above. 

Analytical approaches hke the HPM and ADM that lead to time-power series 
[8-11] are completely useless for a serious study of nonlinear dynamics because 
they fail to provide the overall features of the problem that any physical 
or ecological application requires [1-7]. Other approaches hke the VIM [13] 
also yield expressions that are just valid in a meaningless neighbourhood of 
origin [3]. 

3 Homotopy perturbation method for a prey-predator model 

In a recent paper Rafei et al [12] proposed the application of the HPM to the 
simplest model of prey-predator interaction. After a rather tedious develop- 
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ment of the main equations, their particular implementation of the HPM leads 
to a series solution of the population of the species. It is difficult to justify 
the application of the HPM in this way because one obtains the same solution 
more easily and straightforwardly by direct expansion of the nonlinear differ- 
ential equation in time-power series [9]. Besides, we have recently shown that 
the time-power series developed by Chowdhury et al [11] by means of an alter- 
native implementation of the HPM completely fails to yield the main features 
of the population dynamics [1]. The ADM also leads to a time-series solutions 
although in this case they are slightly different from the correct ones [8]. 

The purpose of this section is to analyze the results of Rafei et al [12] in 
terms of what one expects from an approach designed to solve problems of 
population dynamics. The results of this analysis applies also to the other 
methods already mentioned above [8,9,11]. 

Rafei et al [12] chose the prey-predator system 
dx(t) 



dt 
dyjt) 
dt 



x{t)[a - by{t)], a,b> 
-y{t)[c- dx{t)], c,d>0 (7) 



where x{t) and y{t) are the populations of rabbits and foxes, respectively. This 
nonlinear system ([71) exhibits a saddle point at {xs,ys) = (0,0) and a center 
at {xs,ys) = {c/d,a/b). Besides, the populations obey the following curve in 
the X — y plane: 

In {x^y"") -dx -by = In (xqI/o) - dxo - byo (8) 

where xq and yo are the initial populations at time t = 0. 

The HPM proposed by Rafei et al [12] leads to time-series expansions of the 
form 
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x{t) =xo + xo{a — byo)t + . 

y{t) = yo + yoidxo- c)t + . 



(9) 



that are exactly the same as those obtained earher by Biazar et al [9]. Ob- 
viously, they are suitable about the point (xo,?/o) and will not predict the 
main features of the population dynamics which is what really matters in this 
field [15]. 

Fig. [1] shows the populations for Case I (a = 6 = = 1, c = 0.1, a;o = 14, ?/o = 
18) in a time scale larger than that considered in the earlier studies already 
mentioned above [8,9,12]. We clearly appreciate that all those approaches 
predict a wrong behaviour of the populations. 

Case V {a = b = c = d = 1, xq = 3, yo = 2) is much more interesting. Fig. [2] 
shows the results in the plane x — y. The time series predict a completely 
wrong behaviour and the resulting curve cross itself which can never happen 
as everyone knows [15]. 

Present analysis clearly shows that the HPM [11, 12], the ADM [9] and the 
straightforward time-series expansion [8] are completely useless for a reason- 
able prediction of the evolution of even the simplest prey-predator systems. 

4 Simple epidemic model for non— fatal disease 

In a series of papers Biazar [17] and Rafei et al [18, 19] discussed the applica- 
tion of VAPA to the simplest epidemic model of a non-fatal disease, commonly 
called SIR and attributed to Kermack and McKendrick [20]. Those authors 
applied ADM [17], VIM [19], and HPM [19], and obtained the Taylor expan- 
sions of the number of susceptibles, infectives, and removed about the initial 
time t = 0. In those papers the authors verified that the three methods yield 
exactly the same expansions, and, consequently, they showed the same figures 



6 



for the evolution of the epidemic. 



However, Biazar [17] and Rafei et al [18, 19] did not try to simulate any real 
epidemic situation and did not explain the reason for choosing their particular 
values of the model parameters and initial conditions. 

The purpose of this section is to investigate to which extent the analytical 
expressions derived by those authors are useful for a reasonable prediction of 
the behaviour of an epidemic within the realm of the rather oversimplified SIR 
model. We first answer some relevant questions about the epidemic dynamics 
by means of the exact solution. Then we choose particular model parameters 
and study if the approximate analytical expressions proposed by Biazar [17] 
and Rafei et al [18, 19] are suitable for answering those questions. Finally, we 
give our opinion about the utility of those analytical expressions. 

The nonlinear ordinary differential equations that predict the evolution of the 
epidemic according to the SIR model are: 



^ = —dxy 
dt ^ ^ 

dy n 

— = (3xy - 7?/ 

1^7. (10) 

where x is the number of susceptibles, who do not have the disease but could 
get it, ?/ is the number of infectives, who have the disease and can transmit it 
to others, z is the number of removed, who cannot get the disease or transmit 
it, and /3,7 > are model parameters that determine the epidemic evolution. 

Biazar [17] proposed to solve the SIR equations ( ITOl) by means of the ADM 
and merely obtained the Taylor expansion of the solutions about t = 0: 

oo oo oo 

x{t) = y: x,t\ y{t) = y: y,t\ z{t) = Y z,t^ (11) 

j=0 j=0 j=0 
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where xo = x{0), yo = y{0), zq = z{0) are the initial conditions. Obviously, 
if we substitute the series (ITT!) into the differential equations (fTOj) we easily 
obtain the coefficients Xj, yj and Zj in terms of the initial conditions and the 
model parameters without recourse to any elaborate method like the ADM. 

Later, Rafei et al applied the VIM [f8] and the HPM [19] and obtained the 
same Taylor series. In fact, the three papers show almost the same partial 
sums of the series flTT]) and for that reason display approximately the same 
figures for the discussion of the results [17-19]. 

The SIR model is the simplest description of the evolution of an epidemic of 
a non-fatal disease. It does not consider birth or death of the individuals and 
therefore the total population is constant 

x{t) + y{t) + z{t) =xo + yo + zo (12) 



Besides, the model is so simple that it can be solved exactly and one obtains 



y{x) =yo + Xo- X + - In — 

z{x)=Zo- - In— (13) 
P Xo 

However, if one needs the solutions in terms of t one has to solve the following 
integral numerically: 

1 f dx' 



f^i x' {yo + xo-x' + ^ In ^ 



There are several important questions about the epidemic that one would 
like to answer, for example, if we introduce a small number of infectives in a 
population of susceptibles, will the number of infectives increase, causing an 
epidemic, or will the disease fizzle out?. Assuming there is an epidemic, how 
will it end?, will there still be susceptibles left when it is over?. How long will 
the epidemic last? [21]. Biazar and Rafei et al [17-19] merely showed some 
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figures that are meaningless if tliey did not answer any important question 
about tlie epidemic. 



Tlie exact solution enables one to answer the questions above. As an example 
we analyze how the epidemic dies out. Notice that dx / dt = dy / dt = dz / dt = 
when y = 0, which according to the exact solution takes place when x = xl 
that is a solution of 

yo + Xo - XL + ^In— = (15) 

P Xo 

We may say that the epidemic is over when the number of infectives is similar 
to the one we had at the beginning: yover = Ho] in this case we have 

Xo - Xover + ^ In ^^^^ = (16) 

P Xo 

Biazar and Rafei et al [17-19] chose the following example 

Xo = 20, yo = 15, zo = 10, p = 0.01, 7 = 0.02 (17) 



What we first appreciate here is that the number of infectives is rather large: 
1/3 of the total population. This situation does not define the beginning of 
the epidemic but and advanced stage of it. However, we analyze the results 
anyway. 

The fifth-degree polynomials proposed by Biazar and Rafei et al [17,19] do not 
give acceptable solutions for y = ot y = yo] therefore, they cannot answer 
the relevant questions raised above. In order to find the cause of the failure of 
the time-power series, we show y{x) and z{x) in Fig. [31 We appreciate that 
the time series do not follow the exact curve in the region were we expect to 
find XL = 5.02 x 10"^ and Xover = 0.000908. 

If Xo > 7//? the number of infectives increases up to a maximum at Xm = 7//? 
and then decreases. The figures shown by Biazar and Rafei et al [17-19] are 
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restricted to a rather small time interval and do not show this maximum. 
Besides, they suggest an ever increasing number of infectives and do no reveal 
anything relevant about the future of the dynamics of the epidemic process. 
For the particular model parameters f|T7|) we obtain Xm = 2 and ym = 28.39. 

Although the power-series approach is unable to show the overall picture 
of the epidemic dynamics it does not perform too poorly in the case of the 
conveniently chosen model parameters (fT7|) . One expects that the range of 
utility of the time-power series decreases as the model parameters P and 7 
increases. For example, if we choose 

xo = 20, 2/0 = 4, ^0 = 10, /3 = 7 = 1 (18) 



the situation is considerably worse than in the preceding case. Fig. H] shows 
the exact curves y{x) and z{x) and those predicted by the fifth-degree power- 
series approach. We appreciate that in this case the time series fails to give 
any reasonable account of the epidemic dynamics. 

The results derived above clearly show that the ADM, VIM and HPM imple- 
mented by Biazar [17] and Rafei et al [18,19] are unsuitable for the description 
of the epidemic dynamics. The reason is that the time-power series is a local 
approximation valid in a relatively small neighbourhood of the initial stages 
of the epidemic. Consequently, it cannot provide the long-term behaviour of 
the infectious process that is needed for understanding its future evolution. 
As we have already seen, the analytical expressions provided by Biazar [17] 
and Rafei et al [18, 19] do not allow us to answer the most relevant questions 
about the epidemic dynamics, even in the case of an oversimplified model with 
exact analytical solution. It is most probable that their performance may be 
even poorer in the case of a more elaborated model. 

In conclusion, we do not recommend a health-care system to rely on the ADM, 
VIM or HPM time-power series to cope with an actual epidemic emergency. 
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5 Conclusions 

In our opinion most of tlie recent applications of VAPA are responsible for the 
poorest scientific papers ever published. They are spreading like an epidemic 
overcoming the refereeing mechanism of the journals. I believe that the main 
reason for it is that the authors of such papers referee themselves and accept 
the manuscripts that give them a considerable number of citations. Unfortu- 
nately, it seems that some editors are accomplices of this situation favouring 
such papers and banning criticisms of them. For example, some editors think 
that a prey-predator model that predicts a negative number of rabbits is a 
valuable scientific contribution to the journal [5]. The reader may find some 
more examples elsewhere [1-7]. 
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Fig. 1. Exact (solid, e) and series (dashed, s) populations for Case I 
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Fig. 2. Exact (solid) and approximate (dashed) populations for Case V in the x — y 
plane 
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Fig. 3. Number of infectives y and recovered z in terms of the number of susceptibles 
X for the model parameters (|17p . The sohd and dashed hues correspond to the exact 
and approximate results, respectively. 
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Fig. 4. Number of infectives y and recovered z in terms of the number of susceptibles 
X for the model parameters (jlSp . The solid and dashed lines correspond to the exact 
and approximate results, respectively. 
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